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Abstract 

We present a variational approach to general Lienard type equation in order to linearize it and 
as an example Van der Pol oscillator is discussed. The new equation which is almost linear is factor- 
ized. The point symmetries of the deformed equation is also discussed and the two dimensional Lie 
algebraic generators are obtained. 
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1 Introduction 

Investigation of partial and ordinary nonlinear differential equations is one of the popular recent problems 
in mathematical physics. A typical example of a relaxation oscillator is known as Van der Pol limit cycle 
PP. Van der Pol equation can be adapted to study the Lienard equation. This paper mainly devoted to 
study Lienard type equations whose dynamic properties have been widely studied due to the applications 
in many fields such as physics, mechanics and engineering [H [21 El SI [SI El EJ [5] - This kind of general 
equations and their properties are studied within many methods. Factorization method [91 is one of 
these topics and it can be applied to a class of nonlinear differential equations including commonly 
Lienard type equations (TUl US U2]- Beyond the factorization method, Lie group approach which is an 
intcgrability method is a powerful tool for obtaining solutions and investigating the symmetry properties 
of the differential equations [131 [TS]. In this paper, variational approach and factorization of the 
system are discussed in section 2. In addition, section 3 involves a short discussion of the Lie symmetries 
of the system, the results are concluded in section 4. 

2 Linearizing and factorization 

During the process of deriving Euler's equation one writes 



y(x,a) = y(x,0) + ay x {x). 



(1) 



If we write a difference equation of the form 



A/ 



f[x, y[x, a), y (x, a)} - f[x, y(x, 0),y (x, 0)] 



(2) 




(3) 
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and then continue with introducing the variation of / as 

5 f=% s » + w Sy ' (4) 

where Sy = oty\ and Sy — oty l . Thus the variation of / is expressed in terms of constant a. Thus, this 
operation can be applied to a differential equation called as general Lienard equation: 

x + F(x{t))x + x = (5) 

where the dots correspond derivatives of x(t) with respect to t. Now apply the variational operator to ([5]) 
that is defined by Sx n = nx n ~ 1 5x — y n . Here the differential and variational operators are commuting 
such that 5jfi = 4t5x. In this condition (|5|) becomes 

y + F(x(t))y + S(F(x(t)))x + G{t)y = 0. (6) 

As an example, we can discuss Van der Pol type equation where F(x(t)) = e(l — x(t) 2 ). In this case ([5]) 
will be introduced as a complex differential equation [IB] 

z-e(l-\z\ 2 )z + z = (7) 

where the complex variable z = zi+iz 2 and the equation admits the translation symmetry z(t) —¥ z(t+S), 
phase symmetry z(t) — > e ie z(t), spatial symmetry z(t) — > —z(t), conjugation symmetry z(t) — > z*(t). (J7]) 
can be written in a general form: 

z-e(l- \z\ 2 )z + G(t)z = 0. (8) 
where G{t) is an unknown function. If the variational operation is applied to ([8]) 

5z - e{l - \z\ 2 )Sz + ez 6\z\ 2 G(t) 6z = (9) 



is obtained and © turns into the form 

Y + e(\z(t)\ 2 - l)Y + ei{t)Y 2 + G(t)Y = (10) 

here 5z\ = Y(t) is used, we assume that Sz2 = in the equation given above, according to non-analytical 
zz* . Henceforth, we can discuss the factorization of (fTU)) . We use a homographic transformation D(a, (3, 6) 
that is generally given by 

(y, t) -> (y, T) : y = °® ^ + , aS - h + 0, T = 8(t), p{t) = W) (11) 

and makes the ordinary differential equation (|10p is form invariant. Under this transformation D(a, /3, 9), 
(UUJ) turns into 

d 2 Y I ( a 0" /9 A 

- £ ^ + f^+2 £ ^- £ (l-, 2 )^^ ^ 
z \ a a 10 z 

+ (ff" + ei/3 2 - e(l - z 2 )/?' - G(t)) = 0, afl' 7^ 

a L v z 

In order to transform (fT2]) into an equation below that is addressed as canonic form of type I in |17j : 

^ - 6Y 2 + 6f5 = (13) 

one must write 

/, exa 
6 = — 



(14) 



(Ln a)' = -%{l-x*) + ^) 

j3 = — —(-G(t) + e(l -x 2 )(Ln a)' + (Ln a)" + (Ln a)' 2 ). 
2ex 



In [18] it is pointed out that ()T4|) is essential to obtain (fl3|) and equations in (|14j) must be satisfied. In 
this study we will follow a more general way, there won't be such a restriction given in (|14[) . We will use 



Fi(t) = e(\z(t)\ 2 - 1), F 2 (t) = ez(t) in dTUJ) that can be factorized as 



" MY))(i-fi(Y))Y = (ir.) 



dt J y ' ) \dt 
where f± and f 2 satisfy the relations: 

Fx = -(fi+f2+Y$£) (16) 
F 2 Y + G = hh- 

We want to show that if (fTUj) admits the factorization for /i and f 2 have the form 



fi = ai^/FzY+iVG) (17) 

h = cr x {J¥?y -iVG) 



where a is a constant and can be found as a = ±iy § by comparing (fTT|) and (fTr?|) . Then one can find a 
relation between F\ and G using the first equation in (|16l) 



Fi = i^v 7 ^. (18) 



The second bracket in (fT5l) allows us to write 



^■vf^vf^ (i9 » 

which is a Bernoulli differential equation and we obtain the solutions of the form |19j 

Y = -rt?(fy[^tf)*} \ ^*) = e^'^' (20) 

In order to derive special exact solutions, assume that (<9 t — f\)Y = ^{t), and it yields 

Y-h[Y)Y = (21) 

^-/ 2 (F)V = (22) 



and this choice will lead to have more general solutions. In [TU] the author selected a constant function 
f 2 (Y) but we want to use a non-constant f 2 which is given by the second formula in (|T7|) . Then, (21) 
can be given by 

Y = h(Y)Y + e ^ hiY)dt '. (23) 
Thus, Y can be derived easily as solving this equation 



The solution is introduced as [TH] 



where t dependent function ip is defined by 



(p 3F 2 (t) 36j 

ip V 2a; (t) uj 



Assume that ip is a constant. Thus, a special solution for uj is given by 

v 2 



= ! (dzbi / I (27) 



4 1 J V 6 

where Ci is an integration constant. Now, let us choose F 2 as 

F 2 (t) = -6e^ VS ^ VSdt ' (28) 
then, using (|28l) and (TT51) . the following expression can be obtained 

i=_£ e =F¥/ t (W a -D««' (29) 

and the solution of 2 is introduced by 

z = ce m (30) 



where c and 9 are some constants. In this case, (|2U)) is satisfied when 



c(c 2 - 1) = T5 



±- 



1-c 2 



(31) 
(32) 
(33) 



Thus, Fx, F 2 and G are expressed by 



Fx = t(c 2 -l) (34) 

F 2 = -c9e m 

G = -(1-c 2 ). 
25 v ; 

3 Symmetries of the deformed equation 

The point symmetry of (TTOj) will be discussed here. As it is stated in [11], the differential equations that 
forms a system 

F l (x,u,u {1} , ...,ti (fe) ) = 0, i = l,...,s (35) 
is invariant under the X which is an infinitesimal generator such that 



XF\ Fi=0 = Q. 



(36) 



So, the infinitesimal generator for (|10j) will have a form 



for a second order equation 



V =.f(x,y,y) 



(37) 



(38) 



where the point transformations in the plane are given by 

x m x + a^(x,y) (39) 
y y + ar](x,y) (40) 

Thus we have 

X(y" -f(x,y,y'))\ y „= f = 0. (41) 
The equation above can also be written [13] 



+ (2Vxy -£xx)y + (Vyy - 2 £ xy)y y £,yy Cfx Vfy 



(42) 



+(Vy - 2 ?.t - 3y £ y )f - {rjx + (vy - £, x )y - y € y )f v > = 
In this case, / = —F\y — F 2 y 2 — Gy and here the derivatives are taken with respect to x and y is used 

instead of Y here. We can group the equations as 

0/) 3 : Zyy = (43) 

(y'f: 773,3,-2^ + 2^=0 (44) 

V ■ 2Vxy - Cxx + £F lx + Fit* = (45) 

y 2 : ZF 2x -F 2 (r) y -2Z x )=0 (46) 

y:^G x + 2F 2 r) - G( Vy - 2£ x ) = (47) 

y° : Jfe* + VxFi - rjG = 0. (48) 

yy : 3G£, = (49) 

Thus we can obtain £ from (|4"3")l 

£ = a(z)y + /3(z). (50) 



On the other hand, from (|44l) . rj takes the form 

rj = W(x)y + S(x). (51) 
W{x) = and a(x) =0 can be obtained by using (g7]) and (gHJ)- If % - 2^ in (gSJ) is used in (gTJ) we get 

2F 2 77 + /?(x)^-G = (52) 

and using (j45|) 

/3(a:) = A + Be Flx (53) 
where A, B are constants. We can obtain S(x) using (|48|) as 

S'(x) = de - "*" + Cae"""* (54) 

where Ci, C*2 are constants and a± can be given in terms of the parameters of (jlO[) : 



2 \ 2 

If we use the values of F\, F^ and G and (|47|) we have 



' « ±J ^ . ,55) 



2F 2?7 - Gfo - 2£ a ) = (56) 
where the derivative of G is zero as it is seen from f|34|) . Then, we can obtain 

a+=-i6 + F 1 , a_ = 16 (57) 



Finally we get the value of 9 which is same with the result obtained in the previous section 

. 1-c 2 



(58) 




Figure 1: The skeleton of (fTUl) for some specific parameters. Y = P, Y = Q. M. = {Y, Y, Y} reduced 
manifold. 



which means that the positive sign should be used for 9 in (33). Thus, two dimensional generators are 
obtained as 

x ' = K + ^S < 59 > 

the commutation between these operators is [Xi,X 2 ] ^ and the group is not closed. 

4 Conclusion 

The flexibility of the variational approach is shown by applying the variational operator to the Lienard 
type equation and linearized Van der Pol differential equation is examined in the complex plane. The 
deformed differential equation is factorized and solutions are obtained. Moreover, it is shown that the 
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time dependent coefficients of the linearized equation that are, F± and F%, found as pure imaginary and 
complex for the real values of c. Thus, interesting result is the parameter e is obtained as i. It means 
that the symmetry properties of the equation can be changed under the time reversal symmetry within 
the contributions of F\ , F2 and G such that the original equation (J7J is invariant under spatial reversal 
but (jlOp is not invariant under Y — > — Y and it is not invariant under both spatial and e — i, i — > —i 
parameter and t — > — t reversal. Besides factorization, the symmetry of equation (jlOp is examined within 
point symmetries and two dimensional generators are obtained. The group is not closed, however we 
have concluded the parameter 9 of z(t) is the same as that obtained with both the symmetry equations 
and factorization results. 

The skeleton for (10) is a flat surface such that this equation has two-dimensional generators, as is 
shown in figure 1. 
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